76 research outputs found
A Well-Balanced Central-Upwind Scheme for the Thermal Rotating Shallow Water Equations
We develop a well-balanced central-upwind scheme for rotating shallow water
model with horizontal temperature and/or density gradients---the thermal
rotating shallow water (TRSW). The scheme is designed using the flux
globalization approach: first, the source terms are incorporated into the
fluxes, which results in a hyperbolic system with global fluxes; second, we
apply the Riemann-problem-solver-free central-upwind scheme to the rewritten
system. We ensure that the resulting method is well-balanced by switching off
the numerical diffusion when the computed solution is near (at)
thermo-geostrophic equilibria.
The designed scheme is successfully tested on a series of numerical examples.
Motivated by future applications to large-scale motions in the ocean and
atmosphere, the model is considered on the tangent plane to a rotating planet
both in mid-latitudes and at the Equator. The numerical scheme is shown to be
capable of quite accurately maintaining the equilibrium states in the presence
of nontrivial topography and rotation. Prior to numerical simulations, an
analysis of the TRSW model based on the use of Lagrangian variables is
presented, allowing one to obtain criteria of existence and uniqueness of the
equilibrium state, of the wave-breaking and shock formation, and of instability
development out of given initial conditions. The established criteria are
confirmed in the conducted numerical experiments
Three-layer approximation of two-layer shallow water equations
Two-layer shallow water equations describe flows that consist of two layers of inviscid fluid of different (constant) densities flowing over bottom topography. Unlike the single-layer shallow water system, the two-layer one is only conditionally hyperbolic: the system loses its hyperbolicity because of the momentum exchange terms between the layers and as a result its solutions may develop instabilities. We study a three-layer approximation of the two-layer shallow water equations by introducing an intermediate layer of a small depth. We examine the hyperbolicity range of the three-layer model and demonstrate that while it still may lose hyperbolicity, the three-layer approximation may improve stability properties of the two-layer shallow water system
Local error analysis for approximate solutions of hyperbolic conservation laws
We consider approximate solutions to nonlinear hyperbolic conservation laws. If the exact solution is unavailable, the truncation error may be the only quantitative measure for the quality of the approximation. We propose a new way of estimating the local truncation error, through the use of localized test-functions. In the convex scalar case, they can be converted into L loc ∞ estimates, following the Lip′ convergence theory developed by Tadmor et al. Comparisons between the local truncation error and the L loc ∞ -error show remarkably similar behavior. Numerical results are presented for the convex scalar case, where the theory is valid, as well as for nonconvex scalar examples and the Euler equations of gas dynamics. The local truncation error has proved a reliable smoothness indicator and has been implemented in adaptive algorithms in [Karni, Kurganov and Petrova, J. Comput. Phys. 178 (2002) 323–341].Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/41709/1/10444_2005_Article_7099.pd
HERO (High Energy Ray Observatory) optimization and current status
The High-Energy Ray Observatory (HERO) is a space experiment based on a heavy
ionization calorimeter for direct study of cosmic rays. The effective
geometrical factor of the apparatus varies from 12 to 60 msr for protons
depending on the weight of the calorimeter from 10 to 70 tons. During the
exposure for 5 years this mission will make it possible to measure energy
spectra of all abundant cosmic ray nuclei in the knee region (3 PeV) with
individual resolution of charges with energy resolution better than 30\% and
provide useful information to solve the puzzle of the cosmic ray knee origin.
HERO mission will make it also possible to measure energy spectra of cosmic
rays nuclei for energies 1-1000 TeV with very high precision and energy
resolution (up to 3\% for calorimeter 70 tons) and study the fine structure of
the spectra. The planned experiment launch is no earlier than 2029.Comment: LaTeX,25 pages, 19 figure
A Diffuse-Domain Based Numerical Method for a Chemotaxis-Fluid Model
In this paper, we consider a coupled chemotaxis-fluid system that models
self-organized collective behavior of oxytactic bacteria in a sessile drop.
This model describes the biological chemotaxis phenomenon in the fluid
environment and couples a convective chemotaxis system for the oxygen-consuming
and oxytactic bacteria with the incompressible Navier-Stokes equations subject
to a gravitational force, which is proportional to the relative surplus of the
cell density compared to the water density.
We develop a new positivity preserving and high-resolution method for the
studied chemotaxis-fluid system. Our method is based on the diffuse-domain
approach, which we use to derive a new chemotaxis-fluid diffuse-domain (cf-DD)
model for simulating bioconvection in complex geometries. The drop domain is
imbedded into a larger rectangular domain, and the original boundary is
replaced by a diffuse interface with finite thickness. The original
chemotaxis-fluid system is reformulated on the larger domain with additional
source terms that approximate the boundary conditions on the physical
interface. We show that the cf-DD model converges to the chemotaxis-fluid model
asymptotically as the width of the diffuse interface shrinks to zero. We
numerically solve the resulting cf-DD system by a second-order hybrid
finite-volume finite-difference method and demonstrate the performance of the
proposed approach on a number of numerical experiments that showcase several
interesting chemotactic phenomena in sessile drops of different shapes, where
the bacterial patterns depend on the droplet geometries
Experimental Convergence Rate Study for Three Shock-Capturing Schemes and Development of Highly Accurate Combined Schemes
We study experimental convergence rates of three shock-capturing schemes for
hyperbolic systems of conservation laws: the second-order central-upwind (CU)
scheme, the third-order Rusanov-Burstein-Mirin (RBM), and the fifth-order
alternative weighted essentially non-oscillatory (A-WENO) scheme. We use three
imbedded grids to define the experimental pointwise, integral, and
convergence rates. We apply the studied schemes to the shallow water equations
and conduct their comprehensive numerical convergence study. We verify that
while the studied schemes achieve their formal orders of accuracy on smooth
solutions, after the shock formation, a part of the computed solutions is
affected by shock propagation and both the pointwise and integral convergence
rates reduce there. Moreover, while the convergence rates for the CU
and A-WENO schemes, which rely on nonlinear stabilization mechanisms, reduce to
the first order, the RBM scheme, which utilizes a linear stabilization, is
clearly second-order accurate. Finally, relying on the conducted experimental
convergence rate study, we develop two new combined schemes based on the RBM
and either the CU or A-WENO scheme. The obtained combined schemes can achieve
the same high-order of accuracy as the RBM scheme in the smooth areas while
being non-oscillatory near the shocks.Comment: 33 page
- …